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Abstract 

We consider the motion of rigid bodies in a potential fluid subject to certain nonholonomic con- 
straints and show that it is described by Euler-Poincare-Suslov equations. In the 2-dimensional case, 
when the constraint is reahzed by a blade attached to the body, the system provides a hydrodynamic 
generalization of the Chaplygin sleigh, whose dynamics are studied in detail. Namely, the equations 
of motion are integrated explicitly and the asymptotic behavior of the system is determined. It is 
shown how the presence of the fluid brings new features to such a behavior. 

1 Introduction and outline 

This paper considers the motion of rigid bodies in a potential fluid in the presence of certain kind of 
linear nonholonomic constraints. This is motivated by studying the dynamics of underwater vehicles 
with large fins, which, in the first approximation, impose restrictions on the velocity of their central 
points relative to the fluid. 

The free motion of such vehicles in a 3-dimensional potential fluid can be described by the finite- 
dimensional Kirchhoff equations, in which the action of the fluid is reflected via the tensor of adjoint 
masses that depends solely on the body shape. In the presence of a large fin this tensor becomes almost 
singular: one of its eigenvalues is very large. As was shown in [IS], see also [I], section 6.4, in the limit 
the dynamics of the body with the fin is described by a vakonomic system with the nonholonomic 
constraint that prohibits the instantaneous motion in the direction of the first eigenvalue. 

It should be noted however, that the vakonomic description of several mechanical systems with 
nonholonomic constraints gives rise to rather unexpected behavior (see again |I3j). 



In our paper we follow a more direct approach. Namely, we describe the motion of the body with a 
fin by the Kirchhoff equations with a regular tensor of adjoint masses and impose linear constraint (s) 
on the velocities, which serve to model the action of the fin. The reaction forces arising from the 
constraint are included according to the Lagrange-D'Alembert principle, which is widely accepted to 
be physically meaningful. 

Our model for the constraint and the reaction forces can also be derived following the anisotropic 
friction approach for realizing constraints described in \19\ I13j . see also [5]. In this approach one 
considers the unconstrained system under the influence of a viscous frictional term that only acts in 
the direction perpendicular to the fin and is proportional to a parameter, say x- The equations of 
motion are then obtained by letting x — oo. The obtained system is in general different from the limit 
vakonomic system described above. 

The idea of modeling the action of the fin on the fluid with a nonholonomic constraint using 
Lagrange-D'Alembert's principle also appears in |23j. Although it is regarded as an idealized situation 
that is physically unattainable, it is the first approximation for their analysis. The authors use it to 
study the motion of an underwater projectile with tail fins moving at high speed. However, due to the 
high speed of the projectile, cavitation effects appear and the general interaction of the fluid and the 
body motion is not modeled with Kirchhoff's approach. 

We thus believe that the systems that we present serve as a first approximation for the motion of 
an underwater vehicle with a large (or very effective) fin that moves in a fiuid in the regime where 
Kirchhoff's approach is valid. 

In the case of a 2-dimensional body on a plane and 2-dimensional fiuid, another motivation for 
the constraint appears: the body can interact with the plane via a sharp blade. This setting gives 
a hydrodynamic generalization of the famous nonholonomic Chaplygin sleigh problem considered in 
detail in [61 US]. 

Contents of the paper. In Section 2 we review the preliminaries that are necessary for writing 
down the equations of motion for our family of systems. We briefly recall Kirchhoff's equations for 
a rigid body moving in a potential flow and the Euler-Poincare-Suslov equations for nonholonomic 
systems on Lie groups with left-invariant Lagrangian and constraints. Towards the end of the section 
the equations of motion for underwater bodies subject to a Suslov or Chaplygin sleigh type constraint 
are given. We then discuss the necessary and sufficient conditions for the existence of an invariant 
measure in a simple case. 

The problem of the motion of the Chaplygin sleigh moving in a potential flow (with no circulation) 
is treated in detail in Section 3. The reduced equations are written down explicitly for a general body 
shape and their qualitative behavior is determined. It is shown that these equations are Hamiltonian 
with respect to a given bracket but do not preserve a measure with a smooth density in the generic 
case. 

We then continue to show that in the presence of the fiuid, the sleigh generically evolves from 
one asymptotic circular motion to another in the opposite direction, although the limit circles do not 
coincide. Their radii are given in terms of the components of the total inertia of the fiuid-body system. 
For the purpose of concreteness, the added inertia tensor is computed explicitly for an elliptical sleigh 
whose contact point P with the plane is located at the center of the ellipse O, and where the knife edge 
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is not aligned with its principal axes. This allows also to calculate the components of this tensor when 
P does not coincide with O. 

Then, in the general case, the asymptotic behavior (radius and course direction along the limit 
circle) is fully determined by the position of the center of mass. 

In Section 4 the reduced equations of motion for the hydrodynamic Chaplygin sleigh are integrated 
explicitly for a generic sleigh. The angular velocity is integrated to give a closed expression for the 
angle that determines the orientation of the body. Even though the position of the sleigh cannot be 
obtained in a closed form, the distance between the centers of the limit circles is computed explicitly. 

Finally, in the Conclusions we motivate a further study of the hydrodynamic Chaplygin sleigh in 
presence of circulation and/or point vortices. 

2 Preliminaries 

Rigid body motion in a potential flow 

The motion of a rigid body in a potential fluid in the absence of external forces was first described 
by Kirchhoff in 1890, who derived the reduced equations for the evolution of the body that do not 
incorporate the fluid itself. Its presence is instead encoded in an "added inertia" matrix that depends 
on the body shape. 

Kirchhoff's equations can be understood as the output of a two stage reduction procedure. In 
the first stage one gets rid of the fluid variables by virtue of the "particle relabeling symmetry". 
In the second stage one eliminates the body configuration variables by homogeneity and isotropy of 
space. As such, Kirchhoff's equations are the Lie-Poisson equations on the co-algebra se(3)* where the 
Hamiltonian is the total energy of the fluid-body system. See |2U[|12j for more details on this reduction. 

We briefly recall Kirchhoff's equations in order to introduce the notation needed for the rest of the 
paper. For a derivation of these equations obtained by balancing the momentum of the body with the 
forces and torques exerted by the fluid see Lamb's classic book on hydrodynamics, [16j, which presents 
a thorough discussion of the problem and remains to date a key reference in the subject. 

We adopt Euler's approach to the study of the rigid body dynamics and consider an orthonormal 
body frame that is attached to the body. This frame is related to a fixed space frame via an attitude 
(or rotation) matrix g{t) £ SO{3). 

Let V(t) G M^, u:{t) G be the linear velocity of the origin of the body and its angular velocity, 
both vectors are written with respect to the body frame. We then have 



where the components of x(t) E are the spatial coordinates of the origin of the body frame at time 
t, and 



The configuration of the body is completely determined by the pair {g{t),x{t)), an element of the 



u{t) = g-\t)g{t) 



Y{t)=g-\tMt) 



(2.1) 
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Euclidean group SE{3). In this way ^ := (cj, V) is thought of as an element of the Lie algebra 56(3) 
that is identified with via the bracket 

[{u, V),{r,,V)] = {uxr,,uxV-rjxY), 

where "x" denotes the standard vector product in M'^. 

The kinetic energy of the body is given by (see, e.g., [T7] ) 

1 3 / 3 ^ ^ \ 

^13= 2^ \ ^'"'i + ^'^B^i^j + '^^B^iUj , (2.2) 
i=l \ j=l J 

where m is the total mass of the body and the constants Xg and /Cg, i,j = 1,2,3, depend on its 
shape and the mass distribution. Here the 3x3 matrix Xg is the usual inertia tensor of the body with 
respect to the chosen frame. If the origin of the body frame is at the center of mass, then /Cg = 0. For 
convenience we introduce the 6x6 symmetric matrix 

f Zb K.B 
^^■=\1CI ml 

(I denoting the 3x3 identity matrix) that defines Tb as a quadratic form on se(3). 
Next, the total energy of the fluid is given by 

where p is the (constant) fluid density, u is the Eulerian velocity of the fluid, and the integration takes 
place over the region occupied by the fluid. 

We assume that the fluid motion takes place in the boundless region of that is not occupied by 
the body. We also assume that the flow is potential, with zero circulation, and is solely due to the 
motion of the body. Under these hypothesis it is possible express Tjr as the quadratic form (see jl6j): 

1 ^ 

= - ^ (M'^^ViVj + r^LOiOJj + 2JCjviUJj), (2.3) 
*,i=i 

where j-,X^ and /C^-, i,j = 1, 2, 3, are certain constants that only depend on the body shape. They 
are referred to as added masses and are conveniently written in 6 x 6 matrix form to define the added 
inertia tensor: 

where Xjr, /Cj^, and Mj^ are the corresponding 3x3 matrices. One can show that the matrix Ij- is 
symmetric. 

In the absence of potential forces, the total energy of the fluid-body system is T = Tg + Tj 
and defines the kinetic energy Lagrangian £ : T{SE{3)) — )• M. The motion of the body in space is 
determined by the geodesic motion with respect to the Riemannian metric defined by C. 
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In view of (2.2) and (2.3), we can write the Lagrangian = Tg + Tjr in terms of the hnear and 
angular velocities of the body (written in the body frame) and this expression does not depend on 
the particular position and orientation of the body, i.e. is independent of (5',x). Thus C is invariant 
under the lifted action of left multiplication on 5-^(3). This symmetry corresponds to invariance under 
translations and rotations of the space frame. The reduction of this symmetry defines Euler-Poincare 
equations on the Lie algebra se(3) or, in the Hamiltonian setting, the Lie-Poisson equations on the 
coalgebra se(3)*. The latter are precisely Kirchhoff's equations that are explicitly written below. 

Define the reduced Lagrangian L : se(3) — )■ M by 

m = \eK, 

where ^ = {u:, V) G M'^ x is thought as a column vector and the matrix I = Ig + Ijr. An element /i 
in the co-algebra se(3)* will be represented as a pair /i = (k, p) G M'^ x and will also be thought of 
as a 6 dimensional column vector. Its action on = (o^, V) is defined by 

(/x,0 = k.^ + p.V, (2.4) 

where " • " is the standard Euclidean scalar product. With this identification, the Legendre transform 
defines the mapping between se(3) and se(3)* given by /i = I^. Explicitly we have = (k, p) where 

k = (Xb+Xf)cj + (/Cb + /C^)V, p = mV + TW^V + (/C^ + /C^)^. (2.5) 

In classical hydrodynamics k and p are known as "impulsive pair" and "impulsive force" respectively. 

The reduced Hamiltonian H : se(3)* — )• M is given by H{iji) = and the corresponding 

Lie-Poisson equations /i = adj_i^/i are then (k, p) = ad*^ v)^^' P)- This gives the Kirchhoff equations 

k = kxct; + pxV, p = pxcj. (2.6) 

In the absence of the fluid (p = 0) we have Ijr = 0. Then, choosing the origin of the body axes at 



the center of mass, we obtain /Cg = and k = Xgo;, p = mV. As a consequence, the equations (2.6) 
decouple, and we recover the well known fact about the motion of the body in vacuum: the center of 
mass moves at constant velocity, whereas the body rotates freely according to the Euler equations. It 
is also well known that in the presence of the fluid this is no longer true, that is, the fluid couples the 
translational and rotational modes of the motion. 



equations (2.1 ) 



Given a solution of (2.6), the motion of the body in space is describing by solving the reconstruction 



In the absence of circulation, the description of the motion of the body in the two-dimensional 
case is obtained in an analogous way. The configuration space for the body motion is SE{2) and we 
ultimately get Lie-Poisson equations on se(2)*. This time we write ^ S se(2) as ^ = (w, V) € M x M^, 



and n € se{2)* as = (A;,p) e M X The pairing between /i and ^ is the analog of \2A\ . Then all of 
the above discussion for the three dimensional case remains true by simply inserting the appropriate 
definition of the matrices % and Ijr that relate the column vectors (A:,p) and (cj, V). These are 3x3 
matrices given by 

X + m{a? + })^) —mb ma \ / t r \ 

ls=\ -mb m , 11-^ = ( jvT ) ' (2-7) 

ma m / \ F ^ / 
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where m is the mass of the body, (a, b) are body coordinates of the center of mass, and X is the moment 
of inertia of the body about the center of mass. This time Ijr is a scalar, ICjr is a two dimensional row 
vector, and A4jr is a 2 x 2 matrix. As before, the elements of Ijr depend solely on the body shape. 

The Lie-Poisson equations /i = adj„i^/x are given in components as 

k = V2P1 - V1P2, 
Pi=ujp2, p2 = -ujpi, 



where p = {pi,P2) and {k, p) = I{u}, V) with I = % + Ijr. The reconstruction equations (2.1 ) take the 
form 

(j) = u), vi = xcoscp + ysiiKj), ^2 = — x sin (/> + y cos (2-8) 

where (j) is the rotation angle between the space and the body frame and (x, y) are spatial coordinates 
of the origin of the body axes. 

The Euler— Poincare— Suslov equations. 

We have seen that Kirchhoff equations for a rigid body in a potential fluid are Lie-Poisson equations 
on se(3)* (se(2)* in the two dimensional case) corresponding to a pure kinetic energy left invariant 
Lagrangian. We are interested in adding left invariant nonholonomic constraints to the system. The 
resulting reduced equations, that are consistent with Lagrange-D'Alembert's principle that states that 
the constraint force performs no work during the motion, are the so-called Euler-Poincare-Suslov (EPS) 
equations. We will write these equations explicitly. 

In general, a nonholonomic system on a Lie group G with a left invariant kinetic energy Lagrangian 
and left invariant constraints is termed an LL system. Due to invariance, the dynamics reduce to the 
Lie algebra q, or on its dual q* if working with the momentum formulation. 

The reduced Lagrangian L : g — t- M defines the inertia operator I : g — t- g* by the relation 

where (•, •) denotes the duality pairing. The reduced Haniiltonian, : g* — t- R, is then given by 

HifJ') = ^(At,rV), for fieg*. 

The constraints can be expressed as the annihilator of independent fixed co- vectors G g*. We 
say that an instantaneous velocity ^ G g satisfies the constraints if 

{ui,O=0, i = l,...,n, (2.9) 

where n is the number of constraints. The constraints are nonholonomic if the set of vectors ^ G g 
satisfying the above condition do not span a subalgebra of g. 

The reduced EPS equations on g* are given by, see e.g. jl], 

n 

jjL = adj-i^/i + y^Ait/j, (2-10) 

i=l 

where the multipliers Aj are certain scalars that are uniquely determined by the condition that the 



constraints (2.9) are satisfied. 
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Underwater rigid body with a left-invariant nonholonomic constraint 



We will be interested in the case G = SE{3) with the coalgebra se(3)* = (k, p) and n = 1, which corre- 
sponds to the motion of an underwater rigid body subject to a linear, left invariant and nonholonomic 
constraint 

a-cj + F-V = 0, (2.11) 

a, F being some constant vectors in the body frame. The constraint is nonholonomic provided that the 
set of vectors (a;, V) that satisfy the above condition do not form a subalgebra of se(3). 



Then, in view of (2.6), the EPS equations (2.10) become 



k = kxci; + pxV + Aa, 
p = p X cj + AF, 



(2.12) 



where (k, p) = I{uj, V) as described by ( 2.5 ), the total inertia operator is I = % and the multiplier 



subspace D C se(3) defined by (2.11) or on its image in se(3) 



A is uniquely determined by the constraint (2.11). The system thus describes dynamics on the linear 



When both vectors a, F are nonzero, it is difficult to present a mechanical interpretation of the 
constraint. 

So, we consider two special cases: 

1) F = 0. That is, the constraint is only on the angular velocity: a • a; = 0. Then the equations 



(2.12) represent a hydrodynamic generalization of the classical Suslov problem, see e.g., [2l]. The latter 



describes the motion of a rigid body about a fixed point in presence of this constraint. 

If the origin of the body is at its mass center (/Cg = 0) and in the added masses the translational 
and rotational components are decoupled {ICjr = 0), then, according to (2.5), 



k = (Ib + I^)uj, p = mV + M^Y, 



and the system (2.12) takes the closed form 

k = k X (Xb + Ij^y^k + p X {ml + Mj^y^p + Aa, 
p = p X {Is + Ijr)-^k. 

If we interpret the momentum p as a vector fixed in space, the above system has the same form as 
the equations of the Suslov problem in the quadratic potential (Clebsch-Tisserand) field U = ■ Ap, 
A = {ml +Aij^)^^ , which was studied in detail in [T31[25] under the assumption that a is an eigenvector 
of (Xg+Ijr)"^, see also [llj. The latter condition is important to guarantee the existence of an invariant 
measure (see the discussion below). 

Apparently, for the general tensor + Ijr the system (2.12) with the constraint a • = has not 
been studied before. 

2) a = 0. That is, the constraint is only on the linear velocity of the body, and the system can 
be interpreted as an underwater version of the 3-dimensional Chaplygin sleigh, a rigid body moving in 

under the condition F • V = 0. This mechanical setting is regarded as an approximate model of an 
underwater vehicle with a big fin, as it was mentioned in the Introduction. 



7 



Existence of an invariant measure. It is natural to ask whether the equations (2.12), (2.11) 



possess a smooth invariant measure. This problem has been considered for general EPS equations 
on Lie algebras of compact groups by Kozlov jT^], by Jovanovic fTU] in the non-compact case, and 
by Zenkov and Bloch [2B] for systems with nontrivial shape space. Following |T^, the necessary and 
sufficient condition for equations (2.10), (2.9) to have such a measure in the case n = 1 is that the 
constraint covector u = ui ^ Q* satisfies 



1 



adif- 



c G 



where T G g* is defined by the relation (T, ^) = trace(ad^), ^ G 3. 

Since the group 5-E(3) is unimodular, in our case T = 0, and in case of the generic constraint (2.11 ) 
the above condition becomes 



(a X U + F X W, F X U) = c(a,F), c G M, (2.13) 

where (U,W) :=I-^(a,F). 

This condition is easily analyzed in the special cases that we considered (either F = or a = 0) 
under the assumption that both /Cg = and fCjr = 0. One can then show that (2.13) is equivalent 
to asking that a is an eigenvector of (Xg + Xjr)^^ in the case F = 0, or to the condition that F is an 
eigenvector of ^ = {ml + M.jr)~^ in the case a = 0. 



3 The hydrodynamic planar Chaplygin sleigh 



We now consider in detail the two-dimensional version of (2.12), which corresponds to the coalgebra 
se{2r 



(A;, p) and the nonholonomic constraint 

aoj + Fivi + F2V2 



0. 



Assume that a = and choose, without loss of generality, Fi 
f 2 = and the reduced equations of motion are 



0. Then the constraint takes the form 



k = V2P1 - V1P2, 
Pi = UJP2, P2 = -ujpi + A. 



(3.14) 



I(w, V)-^, with the 3x3 tensor 



0. 



Here the column vectors (A;, p)-^ and (w, V)-^ are related by {k, p)-^ 
I = 1/3 + Ij-. The multiplier A is determined from the condition V2 = 

In the absence of fluid {Ijr = 0) the equations (3.14) become the classical Chaplygin sleigh problem, 
which goes back to 1911, [B], and describes the motion of a planar rigid body with a knife edge (a blade) 
that slides on the plane. The nonholonomic constraint forbids the motion in the direction perpendicular 
to the knife edge. T he asymptotic motions of the sleigh on the plane are straight-line uniform motions, 
see IlEKIg]. 

In the presence of the potential fluid {Ijr ^ 0) the added masses define a more general kinetic 



energy, and the system (3.14) describes a hydrodynamic generalization of the Chaplygin sleigh. We 



shall see that this leads to new features in the asymptotic behavior of the body. 
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The total inertia tensor. Introduce the body reference frame {Ei E2} centered at the contact pomt 
between the knife edge and the plane and choose Ei to be parallel to the blade and E2 orthogonal to 
it. This ensures the above constraint V2 = 0- 



While the expression for Ig with respect to the body frame was given in (2.7) for an arbitrary 
body, the expression for the tensor of adjoint masses Ij- can be given explicitly only for rather simple 
geometries. A simple yet interesting case is an elliptical planar body with the semi-axes A > B > 0. 
Assume first that the origin is at the center of the ellipse, but the coordinate axes El, E2 are not 
aligned with the axes of the ellipse, forming an angle 6 (measured counter-clockwise), as illustrated in 
figure 1 (a). 



Center ol mass 





(a) 



(b) 



Figure 1: The elliptical sleigh. The blade makes an angle 9 with the major axis of the ellipse. 
For this geometry, using the formula for the fluid potential given in [16j, one can show that 



Ijr = p-K 



I (A2-B2)2 
' 4 













52 C0S2 e + ^2 gij^2 Q g-^(2^) 



(3.15) 



V 







sin(26i) A^ cos^ (9 + ^2 sin2 B 



The total inertia tensor, I = + Ij-, of the fluid-body system is then given by 



pn 



\ 



—mb 
ma 



—mb 

pn (S^ cos^ 9 + A'^ sin^ ( 



ma 



sin(26l) 



pTT sin(26') m + pni A^ cos^ 9 + siii^ 9) j 



Notice that in the presence of the fluid, if 7^ n^,n G Z, the coefficient I23 = I32 is non-zero. This 
can never be the case if the sleigh is moving in vacuum as one can see from the expression given for 



IjS in (2.7). The appearance of this non-zero term leads to interesting dynamics that are studied below 



and that, to our knowledge, had not been described before in the literature. 

Now, if the origin of (Ei E2) is not in the center of the ellipse (Figure 1 (b)), then the tensor (3.15) 
takes a more general form with (Ij-)i3, (1^)23 non-zero, which can be calculated explicitly and lead to 
the corresponding modification of the total tensor I. 

In the sequel we assume that the shape of the sleigh is arbitrary convex and that its mass center 
does not necessarily coincide with the origin, which leads to the general total inertia tensor 



I 



J 


-L2 


Li 


-L2 


M 


Z 




Z 


N 
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We keep in mind that we expect to see new phenomena due to the presence of the fluid (when Z 7^ 0). 
The tensor for classical Chaplygin sleigh is recovered by setting Z = 0, J = X+m(a^+6^), M = N = m, 
Li = ma, and L2 = mb. 

Note that, in any case, since the total energy of the motion is positive definite, the tensor I has the 
same property. 



Detailed equations of motion. The constraint written in terms of momenta is V2 = (l ^{k, p)'^)3 
0. Differentiating it and using (3.14), we find the multiplier 



(I-i) 



where 



1 



det(I) 



33 



MN - Z"^ 
ZLi + NL2 
-ZL2 - MLi 



V2P1 - V1P2 

UJP2 

-ujpi 

ZLi + NL2 

JN -L\ 
-L1L2 - JZ 



-ZL2 - MLi 
-L1L2 - JZ 
JM 



^2 



A long but straightforward calculation shows that, by expressing u},vi and V2 in terms oi k,pi,p2, 

constraint V2 = 0, one obt; 

LiUJ + Zvi) — Mvi) , 



substituting into (3.14), and enforcing the constraint V2 = 0, one obtains: 

1 



OJ 



D 
1 

D 



' LiUJ + Zvi) {Jijj — L2V1] 



(3.16) 



where we set D = det(I)(I'i)33 = MJ - L^. Note that D>0 since I and I are positive definite. 



The full motion of the sleigh on the plane is determined by the reconstruction equations (2.8), 
which, in our case with ^2 = 0, reduce to 



X = Vl cos I 



y = Vl sm ( 



(3.17) 



The reduced energy integral has 

H {Juj^ + Mvj- 2L2UJV1) , 

and its level sets are ellipses on the (w vi)-plane. As seen from the equations, the straight line i 
{Liu + Zvi = 0} consists of equilibrium points for the system. 



Hence, if Li and Z do not vanish simultaneously, the trajectories of (3.16) are elliptic arcs that 
form heteroclinic connections between the asymptotically unstable and stable equilibria on i (see Fig. 
2 (a)). The phase portrait is similar to that of the classical Chaplygin sleigh except that the line of 
equilibra is no longer the vi axis if Z 7^ 0. 
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Remark. In fact, the reduced 2-dimensional system (3.16) can be checked to be Hamiltonian with 
respect to the following Poisson bracket of functions of to, 



{Fi,F2} :-- 



D 



{Liuj + Zvi) 



dFi dF2 dFi dF2 
duj dvi dvi duj 



The invariant symplectic leaves consist of the semi-planes separated by the equilibria line £ and the 
zero-dimensional leaves formed by the points on this line. The above bracket can be obtained using 
the construction developed in [H]. 



The case Z = Li = 0. We start by considering the most degenerate case when both Li and Z vanish. 
In particular, for the elliptical sleigh this is the case when = and a = in figure 1 (a). Then all 



the solutions of (3.16) are equilibria. Such a sleigh always performs a uniform circular or straight line 



motion whose parameters depend on the initial condition. In this case, the euclidean measure (or any 



smooth measure for that matter) is trivially preserved by equations (3.16). In fact this is the only case 
in which there is a smooth preserved measure as we now show. 



Proposition 3.1. The reduced equations (3.16) possess a smooth invariant measure if and only if 
Li = Z = 0. 



Proof. The general condition (2.13) for a preserved measure can be specialized to our two dimensional 
problem by putting 

a = 0, F = (0,1,0), U = det(I)-^(-ZL2 - MLi,0,0), W = det(I)-^(-LiL2 - JZ, JM - L^, 0). 
One gets the necessary and sufficient conditions 

MLi + ZL2 = 0, L1L2 + JZ = 0. 



The above can be seen as a linear system of equations for Li and Z with non-zero determinant MJ — 
D > 0. Hence this condition can only hold if Li = Z = 0. □ 



7"2 

^2 



The result of this proposition is to be expected from the qualitative behavior of the system that 
was described above in the case where Li and Z do not vanish simultaneously. 



The case Z ^ 0. In this case the equilibrium points correspond to periodic circular motion of the 
body on the plane, and the contact point of the blade (the origin) goes along circles of the radius 



lim vi I lim a; 

titoo titoo 









z 



(3.18) 



whereas the whole motion is an asymptotic evolution from one circular motion to the other one, but 
in opposite directions, as shown in Fig. 2 (b). (Clearly, when Li = the radius r is zero, and the limit 
motions of the body are just rotations about the fixed origin.) 

The preferred direction of rotation is determined by the following proposition whose proof follows 
from a simple linear stability analysis. 



^ A similar observation for some other generalizations of the Chaplygin sleigh was made in [3]. 
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(a) Reduced phase portrait under 
the assumption Li, Z > 0. The sta- 
ble and unstable equilibra are rep- 
resented by filled and empty dots, 
respectively. 

Figure 2: Reduced phase portrait 



(b) Trajectory of the elliptic sleigh 
on the plane. Asymptotic evolution 
from one circular motion to another 
one in opposite directions. The dot 
on the sleigh surface represents its 
center of mass. 

trajectory of the sleigh in the plane. 



Proposition 3.2. Let the line of equilibra £ = {Liu + Zv\ = 0} be parameterized by vi{s) = 
Lis, uj{s) = —Zs, s G M. The equilibra corresponding to s < are unstable, whereas the equili- 
bra corresponding to s > are stable. 

Proof. We simply perform a linear stability analysis. The matrix associated to the linearization about 
the equilibrium vi = —Lis, co = Zs, is 

{L2Z + MLi)Li {L2Z + MLi)Z 
{JZ + LiL2)Li {JZ + LiL2)Z 

It is seen that this matrix has eigenvalues Ai = (corresponding to the continuum of equilibra along 
the line i = {Liw -1- Zvi = 0}), and A2 = —sE with 

E = JZ^ + 2L1L2Z + MLI = {Z, -Li,0) I {Z, -Li,Of > 0, (3.19) 

since I is positive definite. Thus A2 is positive (negative) if s < (s > 0), corresponding to the unstable 
(stable) direction. □ 



In particular, due to the above proposition and (3.18), the balanced elliptical sleigh depicted in fig 



1 (a), with a > and < ^ < ^, will have a limiting motion in the clockwise direction on a circle 
of radius r = p^(^2_^^2°g^^(2e) as t — )• cxd. Notice that these conclusions on the qualitative asymptotic 
behavior of the system are independent of 6, a feature that is reminiscent of the classical unbalanced 
Chaplygin sleigh (see [19] and the discussion at the end of section 4). However, we shall see (Theorem 
4.2) that the distance between the centers of the circles does depend on b. The conclusions of the 
proposition are also illustrated in figure 2. 
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4 Explicit solution and asymptotic data 

The general solution of the reduced system ( |3.16| ) can be written in the form 

uj{t) = A{ata.nh{At) + a cisech{ At)), 
vi{t) = A{l3tanh{At) + fTC2 sech(At)), 

where the constants 



(4.20) 



DZ DLi r-ZL2 + MLi L1L2 + ZJ 

a = — /3 = -^, ci = VD , C2 = VD , (4.21) 

E = JZ"^ + 2ZL1L2 + MLj, 

only depend on the components of the inertia tensor I. Here a = ±1 corresponding to the two different 
branches of the trajectories on the phase portrait. 



Notice that here the denominator E > 0, as shown in (3.19) and that the arbitrary constant A>0 
is related to the energy H of the system by 

2\E 



The motion on the plane in the general case Z ^ 0. In view of (3.17), the angle (f> is calculated 



by integrating the first expression in (4.20), which yields 

(/>i = y Aaci sech{At) dt = 2crci(arctan(e^*) - vr/4). 



'^2 



J Aatanh{At)dt = a ln(cosh(At)) + (po, 



(4.22) 



(pQ being an integration constant. The angles (pi and 4>2 are an odd and an even function oft respectively. 
One can observe that they have quite different behavior: 

r A + + 7r VD{ZL2 + MLi) 

hm (bi = = ic^TT s ' (4-23) 

t—^±oo 2 2 E 

whereas, as t — t- ibc«, the angle (j)2 asymptotically approaches the linear function l{t) = zizAat + (po — 
a In 2. 

The trajectory of the origin on the plane (x, y) is then described by rather complicated integrals, 
which themselves do not say much about its properties. However, it is natural to calculate the distance 
between the centers of the limit circles. To do this, we shall use 

Proposition 4.1. The centers of the limit circles coincide with the limit positions of the material point 
C, which in the body reference frame (Ei E2) has fixed coordinates {0,f3/a) = (0, —L\jZ). 
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Proof. This can be easily obtained from formula (3.18). Alternatively, in view of (4.20), (4.21), the x- 
and y-velocities of the point C are 



xc = ( vi — —OJ ) cos{4){t)) = Aa 



yc = \ vi-^ijj] sin((/)(t)) = Aa 



Z 
D 



Z 



sech(^t) cos((/)(t)), 
sech{At) sin((/)(t)), 



and both tend to zero as t — )• ±00. Hence C coincides with the centers of the circumferences, since 
otherwise its limit velocity would not be zero. □ 



Setting 



!)i + (j)2 as in (4.22), we obtain the following expressions for the components of the 



vector (Axc, At/c") connecting the centers of the limit circumferences: 

a^/D 



Axe 

^yc 



xc dt 
yc dt 



z 

Z 



—00 
00 



A sech( At) (cos ^\ cos (/»2 — sin (^\ sin (j)2) dt, 
A sech.{At) (sin (pi cos (f)2 + cos (pi sin (j)2) dt. 



Since (pi is an odd and (j)2 is an even function of time, the integrals are reduced to 



Axc 
Aye 



Z 



Z 



—00 
00 



A sech(At) cos (p>i cos (p2 dt, 
A sec]i(At) cos (pi sin (p2 dt 



and, if we set T = At, 



Axc 

^yc 



z 



z 



/oo 
sech(T) cos(2ci(arctan(e^) — 7r/4)) cos(a In(coshT) + (po) dT, 
-00 

sech(T) cos(2ci(arctan(e^) — 7r/4)) sin(a In(coshT) + (po) dT, 



—00 
00 



(4.24) 



ci,Q being specified in (4.21). 

It follows that, like the radii of the limit circumferences, the distance d between their centers does 
not depend on the energy, hut only on the components of the generalized inertia tensor, as expectecQ 

The length scale of this distance is given by the ratio y/D/\Z\ and it depends parametrically on the 
dimensionless quantities a and ci. The explicit dependence of the distance d on the parameters a and 
ci is given by the following 



^ In view of the similarity of the reduced hydrodynamic Chaplygin sleigh and the nonholonomic Suslov problem, the 
distance d can be regarded as an analog of the angle between the axes of the limit permanent rotations of the Suslov rigid 
body in space. 
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Theorem 4.2. The square of the distance d between the limit circumferences is given by 



d^ = {Axcf + (Aye)' 



2ttD 



a 



c\ + 



cosh(a7r) — cos(ci7r) 
sinh(a7r) 



(4.25) 



Proof. Assume without loss of generality that = a In 2 in (4.22) and (4.24). This restriction causes 
a rotation of the vector (Aa^c", Ayc"), but does not affect its length. 

)l3 



First, for simplicity, set ci = 0, that is, (I ^)i3 = 0. Then, in view of (4.22), (j)i = 0, and the 



integrals (4.24) give 



Axc + iAyc = 2- 
= 2 
= 2 
= 2 
= 2 



a^D f°° exp(ialn(e^ +6"^)) 



Z 



e'^ + e-^ 



dT = {z = e^} 



aVD exi){ialn{z + l/z)) 

i2^1 







Z 



dz = {u = arctan z} 







Z 



exp(i a[ln(l + tan^ u) — In(tann)]) du 







Z 



exp ( i a In 



sec^ u 



tanu 

(cos -u)-*" (sin u)-*"^^ 



du 



and, similarly, 



Axc - iAyc = 2 



Z 



(cos u)^" (sin ^i)*" du. 



The last integrals have the form of the Euler 5eta- function (see e.g., [7]) 

B{x,y) = [\''-\l-t)y-^dt = 2 r^\cosuf''-Hsmufy-^du = ^^''^^^^^ 

Jo Jo r 



x + y) ' 



r(-) being the Euler Gamma- function with the properties 



1 



z]r[-+z 



vr 



cos(7r2;) ■ 



r(-z)r(z) 



vr 



z sm vrz 



Then (4.26), (4.27) read 



Axe ± i^yc = 2^^^5(1/2 ^ ia/2, 1/2 ^ ia/2) 



aVD T^{l/2=fia/2) 
~^ r(l =FiQ) 



and, in view of (4.29), the square of the distance is 



d^ =iAxcf + (Ay, 



D [r(i-^f)r(^+if)]' _ I)vr/cos2(7rza/2) 



Z2 r(l + m)r(l-ia) 
D 27rsin(7ria) 27rZ) / sinh(7ra) 



ia(l + cos(7ria)) Z'^a \ 1 + cosh(7ra) y Z'^a \ sinh(7ra) 



ia/ sin(7ria) 
2ttD / cosh(7ra) — 1 



(4.26) 



(4.27) 



(4.28) 



T{z + l) = zT{z), zgC. (4.29) 
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which is real and positive for real a. 



In the general case c\ 7^ 0, under the same changes of variables, the integrals (4.24) yield 



Axc =F i^yc = 2 
= 2 



Z 
aVD 



Z 
aVD 

Z 
a^/D 



cos(2ci {u-TT /4) ) (cos n)^^" (sin n)^^" du 

»7r/2 

cos(2ci(u - 7r/4)) (sin(2u))=^^° du = {w = 2u} 



(2^*°) / cos{ci{w-tt/2)) {sin w)^''^dw 



Z 



(2^^ 



cos 



civr 



cos{ciw) (sin^«)^*" dw 



+ sin 



civr 



sin{ciw) (sinit;)^*" dw 



(4.30) 



The last two integrals in the right hand side can be calculated in terms of the Beta function by 
applying the general formulae 

7rcos(^) 



cos(ciu;) (sinu;)'^ dw 



sin(ciii;) (sinii;)'^ dw 



vB 



2 ' 



vr sm 



'Cl 



u being a complex number with a positive real part. Then, after setting 1/ = 1 it ia, (4.30) gives 



Axc =F ^Ay, 



4=F»"7r 



Z (i±ia)5(l + 2i|^,l + =a^)' 



Therefore, using (4.28), (4.29), we find that the square of the distance is given by: 

D7r2 



d^ = (Axc)' + (Aye)' 



Z2(l + a2) 

Z2(l + a2) 
_D7r2_ 

Z2(1 + q2) 

I?7r2a2 



1 



B{1 + 1 + B{1 + 1 



-Cl— ^a 



r(2 + ia)r(2-ia) 



r (1 + r (1 + =^^3^^) r (1 + r (i + 

a2(l + a2)r(ia)r(-2a) 



(r (^^) r (-^4^)) (r (^) r (-^^)) 



Z2 



sin(7r(^i±^)) sin(7r('^)) 



sin(ia7r) 



47rL> / a \ / sin2 {^^) cosh^ (^) + cos2 {^^) sinh 



Z2 



cf + a2 



2 /'QTT^ 
2 y 



sinh(Q;7r) 



The original references for these formulae are [18] and |27| 
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and the last expression simplifies to (4.25). 



□ 



We mention that the formula of Theorem 4.2 is in perfect correspondence with numerical tests. 



The case Z = 0, Li ^ 0. The condition Z = corresponds to the absence of the fluid, or to the 
case when the blade is parallel to one of two specific perpendicular directions in the body frame with 
the following property: if the solid is set in motion parallel to one of these, without rotation, it will 
continue to move in this manner. For the elliptical sleigh, the two directions are precisely the principal 
axes of the ellipse. 

In this case the system reduces to the classical Chaplygin sleigh, whose motion on the plane was 
described in detail in |19j . Namely, this implies a = and the line of equilibria on the phase plane 
{uj,vi) is the axis uj = 0. The trajectory of the contact point on the plane in this case necessarily has 



a return point and, in view of (3.17), is given by the rather complicated integrals. However, the limit 



behaviors of the sleigh are straight-line uniform motions and, according to (4.23), the angle between 
the limit lines i^ 

^(MLi) 
E ■ 



Acj) = aciTT = (Tvr 



Notice that Ac/) and oj have the same sign (the latter does not change throughout the motion). 
In the absence of the fluid, this expression becomes 



(jvr 



\Jm(X + ma?) 



ma 



which does not depend on b. Figure 3 shows the trajectory of the sleigh in this case for different values 



of (see also [19]). Note that 



2lT 



is the number of "loops" that the sleigh performs in its transition 



between the limit straight line motions. 



Conclusions and further work 

A new series of examples of nonholonomic systems has been presented. These are of interest from the 
point of view of applications in the design of underwater vehicles and mechanisms, since, as we have 
mentioned, the constraint can be interpreted as a simple model for a fln. 

From the mathematical point of view, our systems provide a motivation to study the problem of 
nonholonomic geodesies on the group SE{n) with a general left invariant kinetic energy metric. 

The hydrodynamical version of the Chaplygin sleigh that has been considered, provides a new 
example of a simple, integrable nonholonomic system with an interesting asymptotic behavior. Its 
extensive analysis that we have presented can be of interest in the design of control mechanisms, see 

It should be emphasized that we have only considered the case of zero circulation of the fluid. Our 
preliminary studies show that in the presence of circulation the corresponding equations of motion are 

■^Here A(f> is measured from the line of asymptotic straight line motion in the past to the line of asymptotic straight 
line motion in the future in the trigonometric sense. 
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(a) < A0 < 27r 



(b) 27r < A0 < 47r 



(c) 4tv < A4> < 6tv 



Figure 3: Trajectory of the contact point of tlie sleigh on the plane in the case Z = for different values of A0. 
In all these cases a; > 0, so is an increasing function of time and the sleigh turns counterclockwise. 



no longer of EPS type. However, some of the features of the asymptotic motion remain. In particular, 
for certain initial conditions, one has asymptotic evolution from one circular motion to another, as 
before, but the radii of the limit circles are not the same. 

In this spirit, another interesting problem to consider is to couple the motion of the nonholonomic 
sleigh with point vortices. Such a problem (without nonholonomic constraints) has received interest in 
the last years, see for example [221 El 1^ ■ 

We hope to report with progress on the problems described above in the near future. 
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